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Abstract 

We estimate the mass of the lightest neutral Higgs boson h in the minimal super- 
symmetric extension of the Standard Model with universal soft supersymmetry-breaking 
masses (CMSSM), subject to the available accelerator and astrophysical constraints. For 
m t = 174.3 GeV, we find that 114 GeV < < 127 GeV and a peak in the tan/3 distribu- 
tion ~ 55. We observe two distinct peaks in the distribution of mu values, corresponding to 
two different regions of the CMSSM parameter space. Values of < 119 GeV correspond 
to small values of the gaugino mass and the soft trilinear supersymmetry-breaking pa- 
rameter v4o, lying along coannihilation strips, and most of the allowed parameter sets are 
consistent with a supersymmetric interpretation of the possibly discrepancy in — 2. On 
the other hand, values of > 119 GeV may correspond to much larger values of mu2 and 
Aq, lying in rapid-annihilation funnels. The favoured ranges of m/, vary with mt, the two 
peaks being more clearly separated for m t = 178 GeV and merging for m t = 172.7 GeV. 
If the — 2 constraint is imposed, the mode of the m h distribution is quite stable, being 
~ 117 GeV for all the studied values of m t . 
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1 Introduction 



One of the characteristic predictions of the minimal supersymmetric extension of the Stan- 
dard Model (MSSM) is the mass of the lightest neutral Higgs boson h, which is expected to 
be rrih ^ 150 GeV [1]. This is very consistent with the range rrih < 200 GeV that is favoured 
by global analyses of the available precision electroweak data [2] . Various studies have shown 
that the lightest neutral MSSM Higgs boson is very likely to be discovered at the LHC, and 
possibly at the Fermilab Tevatron collider. It is therefore interesting to attempt to refine 
the MSSM prediction for rrih, and to consider what one would learn from a measurement of 
the h mass [3]. 

We study these questions within the constrained version of the MSSM (CMSSM), in 
which the soft supersymmetry-breaking scalar masses mo and gaugino masses rrii/2 are each 
assumed to be universal at some GUT input scale, as are the trilinear soft supersymmetry- 
breaking parameters A . We impose on the CMSSM the available phenomenological con- 
straints from accelerator experiments, astrophysics and cosmology [4, 5] , treating the super- 
symmetric interpretation of the anomalous magnetic moment of the muon, — 2, as an 
optional constraint, and interpreting the WMAP range for the cold dark matter density [6] 
as an upper bound: Q x h 2 < 0.129. 

We base our study on a statistical sampling of the CMSSM parameter space that is 
uniform in the (mi/ 2 ,m ) plane for 100 GeV < mi/ 2 < 2 TeV, m < 2 TeV, |A /mi/ 2 | < 3, 
2 < tan/3 < 58 and ix > 0, assuming initially that m t = 174.3 GeV [7] and discussing 
later other possible values of m t . We began with a random sample of over 320,000 CMSSM 
points: requiring the lightest supersymmetric particle (LSP) to be a neutralino brought 
the number down to somewhat over 260,000. As seen in Fig. CJa), before imposing the 
various phenomenological constraints we find that rrih 1 is distributed between very low 
values < 110 GeV and an upper limit ~ 128 GeV, with a single peak at ~ 120 GeV. The 
drop-off in the count at low rrih is mainly related to our choice of a uniform measure in the 
CMSSM input parameters: because of the logarithmic dependence of rrih on and m , 
low values of rrih only occur at low values of mi/2, rrio and tan (3. (We recall that rrih evolves 
quickly as rrii/2 is increased at low rrii/2 and more slowly at large m^.) The fall-off at large 
rrih is largely due to our choice of 2 TeV as the upper limit on m^. Extending this upper 
limit would slowly push the peak in Fig. Ufa) to the right, and the count at the peak would 
grow rapidly. Once again, because of the logarithmic dependence of rrih on m^, even a 
modest change in the position of the peak would require increasing the upper limit on rrii/2 

1 We use Fortran code FeynHiggs [8] to calculate rrih- 
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Figure 1: The distribution of the mass of the lightest neutral Higgs boson, m^, in the 
CMSSM (a) before applying the accelerator cuts and the WMAP relic density constraint, 
and (b) after applying these constraints. In the latter case, the red (light) histogram shows 
the points favoured by the optional — 2 constraint. We assume m t = 174.3 GeV in both 
panels. 



substantially. 

We next apply a series of constraints, including the LEP lower limit on the chargino 
mass of 104 GeV, b — > 57 [9] and the limit on £l x h 2 that was discussed above, as well as 
the the direct experimental limit on of ~ 114 GeV [10] 2 . The most severe cut on the 
sample, by far, is that due to the relic density, which for most points exceeds the WMAP 
upper limit. When all cuts are applied our sample is reduced to 3075 points, which are 
plotted in Fig. E^b). Most of the range in irth is still available after imposing the various 
phenomenological constraints, as seen in Fig. ^b). However, we see that the distribution of 
nih within this range exhibits significant structures, with peaks at ~ 121 and 117 GeV, 
and a dip at m h ~ 119 GeV. 

In the rest of this paper, we explain the origins of these features, describe the domains 
of the CMSSM parameter space that populate these peaks in the m/, distribution, and 
discuss the effects of imposing the optional — 2 constraint [12] and varying m t . The 
peaked structures in tan/? and reflect different processes that might reduce the density 
of supersymmetric relics x hito the range allowed by WMAP and other observations 3 : 

2 We recall that, although this lower limit may be relaxed in some variants of the MSSM, its value does 
not change for the CMSSM studied here [11], 

3 Note that although we consider Ao ^ 0, the the stop-coannihilation region [13] is beyond the range we 
scan. 
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either coannihilations with sleptons, most importantly the lighter stau: f\ [14], or rapid 
annihilations: XX y i a t ne heavier neutral Higgs bosons A, H [15], or (exceptionally) rapid 
annihilations: XX the lightest neutral Higgs bosons h [16]. 

The structures in the distribution imply that, once m t is better known from Tevatron 
and/or LHC measurements and assuming the CMSSM framework 4 , a measurement of rrih 
at the LHC or Tevatron might enable one to estimate ranges for the values of mx/2,Ao 
and tan f3, even if sparticles themselves are not yet discovered. If sparticles are discovered, 
confronting their masses with the ranges inferred from will be a crucial test of the CMSSM 
framework. 

2 Effects of Phenomenological Constraints on the CMSSM 
Parameter Space 

As already mentioned, we have sampled uniformly the {mi/2, mo) plane for mi/2, mo < 2 TeV, 
|A)/ m i/2| < 3, 2 < tan/3 < 58 and /i > 5 , assuming m t = 174.3 GeV as our default 6 . As 
tan j3 is increased, there is an increasing fraction of sample points that do not yield consistent 
electroweak vacua. Nevertheless, the consistent solutions are distributed quite smoothly in 
tan (3 before applying the accelerator and cosmological cuts, as seen in panel (a) of Fig. |2 
However, after applying the cuts, the distribution in tan/3 is far from uniform, as seen in 
panel (b) of Fig. El The distribution of allowed models is sharply peaked towards large 
tan/3, with a relatively small tail surviving below tan/3 ~ 20. This observation holds for 
both the general sample and the (g M — 2)-friendly subsample 7 , shown as the light (red) 
shaded histogram in panel (b). 

The preference for large tan/3 is an understandable consequence of the interplay of the 
various accelerator and cosmological constraints. For example, the cosmological upper limit 
on the supersymmetric relic density in the coannihilation region imposes an upper limit on 
mi/2 that is significantly relaxed at large tan/3, in particular by rapid XX —>■ A,H annihila- 
tion. Moreover, the funnels due to the rapid-annihilation processes XX ~~ * H,A are broader 
than the coannihilation strips that define the acceptable cosmological regions at lower tan /3. 
We also note that the range of small values of mi/2 that is excluded by the experimental 
lower limit on diminishes as tan (3 increases, and recall that the predominance of high 

4 We also assume that theoretical errors in the CMSSM calculation of m,h can be reduced along with the 
experimental error. 

5 This sign of \i is suggested by even a loose interpretation of g M — 2. 
6 We discuss later the effect of varying mt in our analysis. 
7 We assume here — 2 range from 6.8 to 43.6 x 10 -10 [12]. 
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Figure 2: (a) The distribution in tan/? of the sample points shown in Fig. QJfaj, before the 
accelerator and WMAP constraints are applied, and (b) their distribution after applying these 
phenomenological constraints. In the latter, we see that the distribution of g^ — 2-friendly 
points (coloured red/grey) is similar to that of the total sample. 

tan/3 in satisfying constraints was clearly seen in a likelihood analysis [17] when comparing 
regions of high likelihood for tan (3 = 10 and 50. 

Panel (a) of Fig. |H] displays the allowed points in the (A ,tan(3) plane. We see that 
they gather into three clusters: one centered around A = that extends to small values 
of tan/?, and two at large values of \A Q \ that are concentrated at larger tan (3, particularly 
for A Q < 0. As seen in panel (b) of Fig. these accumulations populate different regions 
of m h . Specifically, the points with m h G (114, 119) GeV, which populate the low-mass 
peak in Fig. ^b), have relatively low values of Aq, most of which are negative. On the 
other hand, points with A < —2 TeV generally have rrih G (119, 122) GeV and points with 
A > 1 TeV have G (122, 127) GeV. Between these wings, there are addditionally some 
low-|v4 | points with G (119, 124) GeV. Thus, the higher-mass peak in Fig. ^b) receives 
contributions from all regions of A . We also see in Fig.^b) that essentially all the low-mass 
points are (g^ — 2)-friendly (shaded red/grey), that only some of the high-mass points with 
Aq > —2 TeV are (g M — 2)-friendly, and that none of the points with A < —2 TeV are 
(g^ - 2) -friendly. 
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Figure 3: (a) The distribution in the (A , tan (3) plane of the sample points shown in Fig. W(b) 
that survive the accelerator and WMAP constraints, and (b) a scatter plot of these points in 
the (A Q ,rrih) plane. The (g^ — 2)-friendly points are coloured red (grey). 



3 Interpretation of Features in the rrih Distribution 

The origins of many of these features can be understood qualitatively by referring to the 
various (m 1 / 2 ,m ) planes displayed in Fig. 2 of [18] for different values of tan/3 and A Q . 
Updated planes for the case tan j3 = 55, whose importance can be seen from panel (b) of Fig. EJ 
and panel (a) of Fig. EH are shown in Fig.|H When tan (3 < 45 and \i > as assumed here, the 
regions allowed by WMAP and the other constraints are essentially narrow coannihilation 
strips that decrease in width as mi/2 increases, terminating when mi/2 ~ 1000 GeV. In most 
of these regions, \A \ < 1000 GeV also, so these points populate the central A ~ region 
in Fig. Efb). Therefore, they provide the majority of the models in the low-mass peak in 
Fig.^b), but also a tail extending under the higher-mass peak, as seen in panel (b) of Fig. El 
These coannihilation strips are also the dominant features for tan/3 = 55 when Aq > 0, as 
seen in panel (a) of Fig. EJby the two examples for A = +mi/ 2 , +2mi/2- 

However, a second class of features is visible in Fig. 2 of [18] when tan f3 > 50, namely 
rapid- annihilation funnels at large m 1 / 2 , as updated in panels (a, b) of Fig. |U for A = 
0, —mi/2, —2m 1 /2- These funnel regions populate the high-mass peakin Fig.^b). The funnels 
are typically broader than the coannihilation strips, and therefore have a larger weighting in 
the const ant- density sampling of the (7^1/2, ^o) plane that we have made in this paper. As 
a result, the larger values of tan/3 have a strong weight in the sample of models surviving 
the accelerator and WMAP constraints that we showed in Fig. El We recall that we retain 
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Figure 4: Regions in the (777.1/2, mo) planes for tan (3 = 55, as calculated for m t = 174.3 GeV 
using the latest version of the SSARD code [19]. Note the narrow coannihilation strips for 
Aq > in panel (a) and the broader rapid- annihilation funnels for A < shown in panel 
(b). 



in our analysis points whose relic density Q x h 2 falls below the range favoured by WMAP [6], 
which typically have slightly lower values of m than along the coannihilation strips, while 
remaining within the region where the LSP is the lightest neutralino \i an d lie inside the 
rapid-annihilation funnels. Restricting our plots to points with Q x within the WMAP range 
would reduce the statistics in our plots, but not alter their basic features. The weight of 
the rapid- annihilation points could be diminished if one used a different sampling procedure, 
e.g., if one gave less weight to regions of parameter space with large and/or m , and 
hence \Aq\, as might be motivated by fine-tuning considerations. However, the 'twin-peak' 
structure of the rrih distribution would survive any smooth reweighting of parameter space. 

The rapid-annihilation funnels are responsible for the dense cluster of models at large 
tan/3 and Aq < —2 TeV in panel (a) of Fig. El which have > 119 GeV as seen in panel 
(b) of Fig. EJ and hence populate the higher peak in the Higgs mass distribution in Fig. [T] 
It is also clear from Figs. [T] and El that the basic feature of a doubly-peaked Higgs mass 
distribution linked to different ranges of A would also survive any smooth reweighting of 
the parameter space. 

As discussed in [18] and seen in Fig. 01 the locations of the rapid-annihilation funnels 
are very sensitive to A , reflecting the sensitivity of rriA,H to this parameter (among others). 
Starting from the A = case where the funnel extends above m^j ~ 1000 GeV for m t = 
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174.3 GeV, the funnel moves to smaller mi/2 as Aq decreases and merges progressively with 
the coannihilation strip. On the other hand, no funnels are visible for A$ sufficiently > 0. The 
WMAP regions for Aq = ±2mi/2 provide points in Fig. El with extreme positive and negative 
values of A , respectively. The different breadths of these regions explain the asymmetry in 
panel (b) of Fig. El in particular. Indeed, for A > these points are simply continuations 
of the coannihilation strips. As seen in Fig. El some of these points are (g^ — 2)-friendly, and 
provide the tail under the high-mass peak in Fig. ^ 

We now consider the (mi/2, A ) planes shown in Fig. El where panel (a) shows the combi- 
nation of all vales of tan f3, and panel (b) shows only models with tan f3 > 50. These panels 
update analogous plots in [18], and display significant differences due to the reduction in m t 
from 178 GeV to 174.3 GeV and improvements in the treatment of vacuum stability require- 
ments. Previously, we had seen a clear separation between 'fins' at A ~ ±1.5rai/2 and a 
'torso' at Aq ~ 0, which has vanished apart (possibly) from a vestigial fin at Aq > 2 TeV that 
is more visible at large tan (3. We also note the appearance of a small 'head' with a 'tooth' 
at mj/2 ^ 150 GeV and A ~ 0, which is due to points with m x ~ m/j/2, whose relic density 
falls within the WMAP range thanks to rapid annihilation through the light CMSSM Higgs 
pole [16]. These points have m x ± very close to the LEP lower limit, and might be accessible 
to the Tevatron. 

We see in panel (a) of Fig. El that the points with < 119 GeV (darker blue/black and 
red/grey colours) cluster at rrii/ 2 ,A < 1 TeV and A > —3 TeV. Almost all these points 
make a supersymmetric contribution to — 2 that could explain the possible discrepancy 
between experiment and the Standard Model calculation based on e + e~ data (indicated in 
red/grey). On the other hand, only a small fraction of the models with nth > 119 GeV (pale 
colours) are compatible with this supersymmetric interpretation of — 2 (pink/light grey). 
As seen in panel (b) of Fig. El all the parameter sets with tan/3 > 50 have rrih > 119 GeV. 
The [g^ — 2)-friendly points are concentrated at mi/ 2 ^ 1 TeV. 

This analysis can be used as a diagnostic tool when the Higgs boson is discovered at 
the Tevatron or the LHC, at least within the CMSSM framework and assuming that m t = 
174.3 GeV. This framework would be invalidated if > 127 GeV. On the other hand, 
if the Higgs boson is discovered with a mass < 119 GeV, one can infer from Fig. Et a ) 
that both mi/2 and Aq must be small, and that supersymmetry is likely to lie along a 
coannihilation strip. On the other hand, if > 119 GeV, supersymmetry may well have 
chosen a rapid-annihilation funnel. 
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Figure 5: (a) The distribution in the (mi/2, A)) plane of the points shown in Fig. U](b) that 
satisfy the accelerator and WMAP constraints, separated into the low-rrih region (blue points) 
and the high-nih region (red points), (b) The distribution in the (m!/ 2 ,^4o) plane of points 
with tan (3 > 50. 



4 Potential Impact of — 2 

We now comment further on the potential impact of imposing the — 2 constraint [12], 
which we treat as optional. We see in Fig. H^b) that this constraint would suppress the 
high-m^ peak, while retaining most of the low-m^ models. The suppression of the high-m/j 
peak is a consequence of the removal of points with large and/or m that would make 
a very small contribution to g^ — 2, many of which are in the rapid-annihilation funnels. A 
similar effect reduces also the upper part of the low-m^ peak, but the coannihilation strips 
would be less affected by the g^ — 2 constraint. On the other hand, as seen in Fig. [21 imposing 
the g^ — 2 constraint would not alter the statistical preference for large tan (3. As we see 
in Fig. 01 imposing the g^ — 1 constraint would disfavour models with large negative Aq, as 
well as many with large positive A , but some models with large tan/? and a small A would 
survive. 



5 Dependence on rrtt 

In all the above, we have assumed that m t = 174.3 GeV [7]. The central value was formerly 
178 GeV [20], and the current central value is m t = 172.7±2.9 GeV [21], following significant 
evolution during recent months. In view of this and the remaining experimental uncertainty, 
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Figure 6: As for panel (b) of Fig. QJ, but assuming (a) m t = 178 GeV and (b) m t = 
172.7 GeV. 

we have also considered the dependence of the above analysis on m t . We recall that Amj ~ 
1 GeV for Am t = 1 GeV in theoretical calculations, and that the parameter regions allowed 
by WMAP vary quite considerably with m t , particularly in the rapid- annihilation funnel 
region, as seen in Fig. 1 of [18]. Specifically, this region moves to smaller (larger) mi/2 for 
smaller (larger) m t . As was already mentioned, the coannihilation strips mainly populate the 
low-mass peak in Fig. Q whereas the high-mass peak is largely due to the rapid-annihilation 
funnels. Accordingly, we would expect these peaks to be more separated at large m t than 
at smaller values. Precisely this effect is seen in the two panels of Fig. |H1 We see in panel 
(a) that the upper peak in Fig. [TJb) shifts upwards by ~ 4 GeV if m t = 178 GeV [20], and 
is very clearly separated from the low-mass peak. Correspondingly, the upper limit on 
increases to 130 GeV. On the other hand, we see in panel (b) of Fig. El that the two peaks 
merge for m t = 172.7 GeV [21], and the upper limit on rrih decreases to 126 GeV. Likewise, 
many of the other features discussed previously in this paper become more (less) pronounced 
for larger (smaller) m t . 

By the time the Higgs boson is discovered, we expect m t to be known with considerably 
greater precision than the present uncertainty 5m t ~ 2.9 GeV. Once m t is known with an 
accuracy 5m t < 2 GeV, and assuming that the accuracy of theoretical calculations in the 
CMSSM can be brought to the same level, there will be no theoretical ambiguity between the 
Higgs mass peaks, and diagnosis of the supersymmetric parameters will indeed be possible 
along the lines discussed in the previous Section. 
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As in the case m t = 174.3 GeV shown in panel (a) of Fig. ^ for m t = 178 GeV the effect 
of imposing the — 2 constraint would also be to remove the high-mass peak, leaving a 
plateau extending from m t ~ 118 GeV to ~ 127 GeV. The low-mass-peak would also be 
reduced, but most of the intermediate plateau for m t = 178 GeV would survive the g M — 2 
constraint. In the case of m t = 172.7 GeV shown in panel (b), there is a more pronounced 
peak at m t ~ 117 GeV and a tail extending to ~ 124 GeV. It is striking that, whatever the 
value of m t , the mode of the distribution is relatively stable at ~ 117 GeV and that the 
upper limit on also remains relatively stable around 125 GeV, if the g^ — 1 constraint is 
imposed. 

6 Conclusions 

We have shown that the available experimental and cosmological constraints on the CMSSM 
allow only a limited range of m^. If m t = 174.3 GeV, this is < 127 GeV without the — 2 
constraint and < 124 GeV if it is imposed. If g M — 2 is not imposed, we find twin peaks in 
the ruh distribution at rrih ~ 117, 121 GeV. The upper bound and the lower peak are quite 
insensitive to variations in m t , whereas the upper peak is sensitive, and merges with the 
lower peak for low m t . Large values of tan (3 ~ 55 are favoured by our analysis, whether the 
— 1 constraint is applied, or not. 

We have also shown in this paper that the mass of the lightest CMSSM Higgs boson may 
be a useful diagnostic tool for identifying the most likely regions of the CMSSM parameter 
space, even if sparticles are not (yet) discovered. This is because the CMSSM is divided into 
distinct coannihilation and rapid-annihilation regions, and measuring could provide us 
with a hint which alternative is chosen by Nature. 
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